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Quantum networks are composed of quantum nodes that 
interact coherently by way of quantum channels and 
open a broad frontier of scientific opportunities*. For 
example, a quantum network can serve as a 'web' for 
connecting quantum processors for computation^^ and 
communication'*, as well as a 'simulator' for enabUng in- 
vestigations of quantum critical phenomena arising from 
interactions among the nodes mediated by the channels^ 
The physical realization of quantum networks generi- 
cally requires dynamical systems capable of generating 
and storing entangled states among multiple quantum 
memories, and of efficiently transferring stored entan- 
glement into quantum channels for distribution across 
the network. While such capabilities have been demon- 
strated for diverse bipartite systems (i.e., N — 2 quan- 
tum systems)^ '^, entangled states with A^ > 2 have 
heretofore not been achieved for quantum interconnects 
that coherently 'clock' multipartite entanglement stored 
in quantum memories to quantum channels. Here, we 
demonstrate high-fidelity measurement-induced entangle- 
ment stored in four atomic memories; user-controlled, co- 
herent transfer of atomic entanglement to four photonic 
quantum channels; and the characterization of the full 
quadripartite entanglement by way of quantum uncer- 
tainty relations*^ *^. Our work thereby provides an im- 
portant tool for the distribution of multipartite entangle- 
ment across quantum networks. Moreover, we show how 
our entanglement verification method, originally devel- 
oped for infinite-dimensional bosonic systems, also applies 
to finite-dimensional quantum spin- 1/2 systems, where 
our results may find application for investigations of en- 
tanglement order in condensed matter systems at ther- 
mal equilibrium** *^. With regard to quantum measure- 
ments, the multipartite entangled state stored in the quan- 
tum memories can be appUed for sensing an atomic phase 
shift beyond the Umit for any unentangled state. 

Diverse applications in quantum information science re- 
quire coherent control of the generation, storage, and trans- 
fer of entanglement among spatially separated physical 
systems' ^. Despite its inherently multipartite nature, entan- 
glement has been studied primarily for bipartite systemPSE^ 
where remarkable progress has been made in harnessing 
physical processes to generate 'push-button' and 'heralded' 
entanglemenPESl as well as to map entangled states to and 
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from atoms, photons, and phonons^llll in this endeavor, 
well-established methods for characterizing bipartite entan- 
glement of discrete and continuous quantum variables have 
been essentiaPSMU 

For multipartite systems, the 'size' of a physical state 
characterized by its density matrix p^ grows exponentially 
with the number of subsystems A^ and makes the entangled 
states exceedingly difficult to represent with classical infor- 
mation. Importantly, this complexity for pat increases the 
potential utility of multipartite entanglement in quantum in- 
formation science, including for quantum algorithms'^ and 
simulatiorPS. Redundant encoding of quantum information 
into multipartite entangled states enables quantum error cor- 
rection and fault tolerant quantum computation'^. The intri- 
cate long-range quantum correlations of quantum many-body 
systems are intertwined with multipartite entanglement in a 
fashion that is a subject of active investigation"''^. In ad- 
dition, mobilizing multipartite entanglement across quantum 
networks could lead to novel quantum phase transitions for 
the network, as for the percolation of entanglement^ . 

Counterposed to these opportunities, the complex structure 
of multipartite entanglement presents serious challenges both 
for its formal characterization and physical realization. Theo- 
retically, there is no unique measure of M-partite entangle- 
ment for M > 2 components of the A^-partite p^v, nor is 
there a complete classification of the types of entanglemenP^l 
Apart from the technical difficulty of preparing pM for a large 
system, there is only a nascent understanding of robust means 
for entanglement verification and characterization that incor- 
porate experimental imperfections and finite measurement re- 
sources. Indeed, there are relatively few examples of labora- 
tory systems that have successfully generated multipartite en- 
tanglement. Most works have considered the entanglement for 
spin systems, notably trapped ions'^ '^, which are applicable 
to the matter nodes of a quantum network. But the method- 
ologies employed for entanglement characterization are prob- 
lematic for the infinite dimensional bosonic systems of the 
quantum channels (e.g., quadrature'^ and number-state''* en- 
tanglement for optical modes). A relevant advance for pho- 
tonic systems has been a posteriori entanglement generated 
from parametric downconversion'^. However, in this case the 
post-dieted state ctm is obtained by a destructive local filter, is 
only a small component of a larger physical state p^v, and is 
not available for subsequent utilization". 

In addition to the generation and characterization of mul- 
tipartite entanglement, an important capability for quantum 
networks is the development of quantum interfaces capable 
of generating, storing, and dynamically allocating the entan- 
glement of matter nodes into photonic (or phononic) channels. 
As illustrated in Fig. [Th, we introduce here such a quantum in- 
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FIG. 1: Overview of the experiment, a, Quantum memories and interfaces for multipartite quantum networks. Entangled state p\^ for 
four atomic ensembles e — {a, b, c, d} is heralded by a photoelectric detection event at detector Dh derived from quantum interference of 
four fields 71 = {ai, fei, ci, di}. After a user-defined delay r, read lasers are applied to the individual ensembles to coherently transform 
the atomic entangled state p\y into quadripartite entangled beams of light p\^ , with each beam propagating through quantum channels 
72 = {0.2, &2, C2, ^2}. Intensity modulators IM„iiic,,cad control the intensities of the writing and reading beams. Inset i, a fluorescence image 
of the laser-cooled atomic samples {a, b, c, d} (Appendix), b, Entanglement generation. A weak write laser is split into four components to 
excite atomic ensembles {a, b, c, d} via parametric interactions C/write. A single Raman scattered photon for four fields 71 is emitted by the 
ensembles and detected by a projective measurement Ilh at detector Dh, which signals the creation of an entangled state pj^,. for ensembles 
{a, b, c, d} (Eq.fTll. c, Quantum-state exchange and entanglement verification. A strong read pulse is sent into the four atomic ensembles, and 
the atomic state p\y is mapped to an entangled state of light p^ for four fields 72 (Eq. [A via quantum-state transfers, (Tread. This entangled 
field state then propagates to the entanglement verification ports, (i) Upper panel for yc-measurement — The quantum statistics {qij^k,i} of 
p\^ are measured with projectors {n^^ } for the modes at detectors Da.b.c,d- (ii) Lower panel for A-measurement — By rerouting the relevant 
fibre-optical connections, we access mutual coherences for p\^ with projectors {11^} from detection statistics Pijki at Da,b.c,d- States \g), 
\s} are hyperfine ground states F = 4, F — 3 of 6S1/2 in atomic Caesium; state |e) is the hyperfine level _F' = 4 of the electronic excited 
state 6P3/2. 



terface for quadripartite entangled states based upon coherent, 
collective emission from matter to light. We present a system- 
atic study of the generation and storage of novel quadripartite 
entangled states of spin-waves in a set of four atomic memo- 
ries, which could be located within one node or distributed at 
multiple nodes, as well as of the coherent transfer of the en- 
tangled components of the material state into individual pho- 
tonic channels. We measure transitions from A/ = 4 to 3 to 
2-partite entanglement and to a fully separable state as a func- 
tion of controlled spin-wave statistics of the atomic memo- 
ries. We also explore the temporal decoherence of the atomic 
entanglement and observe a dynamic evolution of multipar- 
tite entanglement into various subsets of the full quadripartite 
states in a dissipative environment, from fully quadripartite 
entangled to unentangled. 

As illustrated in Fig. [Tl our experiment proceeds in four 

(A) 

Steps (see Appendix). First, in step (1) an entangled state p\^ 
of four atomic ensembles is generated by quantum interfer- 
ence in a quantum measurement"^^ (Fig. flp). Given a photo- 
electric detection event at Dh, the conditional atomic state is 
ideally a quadripartite entangled state p\y = \W)j^ {W\ with 

\W)^ = l[i\sa,9b,9c,9d)+(i'^'\9a.sb,gc,9d))+ 

(1) 
e"^' (l5a> 56, Sc, 9d) + e*"*" Iffa, 9b, 9c, Sd))], 



excitations \Sf) = 



whose single quantum spin- wave \se) is coherently shared 
among four ensembles e = {a, b, c, d}. These entangled 
states are known as ly-states, whose component states are the 
atomic ground states \g^) — {9 ■ ■ ■ g)e and single collective 

'W= T.f=i^ \g---Si--- 9)e, where NA,e 



is the number of atoms in ensemble e. 

After the heralding event at Dh, step (2) consists of stor- 

(A) 

age of pyy in the ensembles without optical illumination for 
a user-controlled time r. At the end of this interval, step 
(3) is initiated with read beams applied to the individual en- 
sembles to coherently transfer the respective entangled com- 



(A) 

ponents of p\^. into a quadripartite entangled state of light 
Pw = 1^)7(^1 '^ia cooperative emissions'* (Fig. lie), where 

\W}j = ^[(|1000)+e*'^i|0100))+e'''^H|0010)+e*^3|0001))]. 

(2) 
This photonic state is a mode-entangled VK-state'"* '^, anal- 
ogous to Eq. [T] which shares a single delocalized photon 
among four optical modes 72 = {021 &2, C2, ^2}. Because 
the mappings of quantum states from the collective atomic 
modes of the ensembles to individual field modes are local 
operations, the presence of quadripartite atomic entanglement 
for p\y (t) can be unambiguously determined by accessing 



the entanglement degree of the photonic state ply {t) via the 
state transfers f/read at time r. 

Finally, in step (4) we verify and characterize the heralded 
entanglement for pyl/ from cor nplem entary measurements of 
photon statistics and coherence'"* '^ (Fig. lie). In particular, 
we consider the reduced density matrix p^ — PoPo + PiPi + 
P>2P>2 containing up to one photon per mode, which leads 
to a lower bound for the entanglement of the actual physical 
State py^ in an infinite dimensional bosonic space. In turn, the 
truncation provides a lower bound inference for the atomic en- 
tanglement in p\^ ■ Here, {po,pi,p>2} are the probabilities 
for the and 1-photon subspaces po.i and the higher-order 
subspace p>2, respectively. As illustrated in the upper panel 
of Fig. fTt, these quantities are evaluated from the probabili- 
ties Qijki for i, j, k, I E {0, 1} photons to occupy the respec- 
tive optical modes 72 = {a2,b2,C2,d2} at the output faces 
of the ensembles via photoelectric detections {H^* } at de- 
tectors Da,b,c,d- The photon probabiUties {po,Pi,P>2} are 
combined into a normalized measure of the degree of statis- 
tical contamination from multiple excitations (0>2 for p\^', 
namely j/c = I ^-y" , with j/c = for a single excitation 

(i.e., p>2 = 0) to j/c = 1 for balanced coherent states'^. Op- 
erationally, we control tJc (thereby, the spin-wave statistics for 

p\y and the photon statistics for p\y) by way of the intensity 
of the write beam. 

We must also quantify the mutual coherences dap 

for the optical modes a,^ e {02,^2,22,^2} of p\^ . 
This is accomplished by measuring the photon probabili- 
ties {piooo,Poioo,Pooio,Poooi} at the outputs of the ver- 
ification (v) interferometer shown in the lower panel of 
Fig. fTt, from which we determine the sum uncertainty 



(n|'^')^) for the collective variables 



{Il[^^} = {m),{W,\}. Here, {\W,).,} is a set of four or- 
thonormal W-states, with phases {/3i, /?2, /^s}?) for the set se- 
lected by way of the actively stabilized paths in the verifica- 
tion interferometer For appropriate choices of {/3i, /32, /^aji;, 

' (c) 

photodetections at Di act as projective measurements 11^ 
of the input pr onto |tVi)^. From our measurement of A, 
we deduce an effective visibility 14ff = Ad for the interfer- 
ence of any two of the four modes {02,62,02,^2}, where 
d — a\g{dap) with < d < j. Hence, for the ideal Ha- 
state in Eq. [2] we have dap — 1/4, so that A = and 
Ves = 1, with associated photon probabilities piooo = 1 and 
Poioo = Pooio = Poooi = 0. By comparison, we obtain 
Piooo ^ 0.97±0.01,poioo ^Pooio ^Poooi ^ O.OlOi 0.003 
in the experiment for ^c ~ 0.04 ± 0.01, as shown in the 
bar plots of the lower panel of Fig. [T];. In contrast, a 
mixed state analogous to Eq. |2]but with no phase coherence 
dap — would result in balanced photon probabilities (i.e., 
Pwoo = Poioo = Pooio = Poooi = 1/4) at the outputs of the 
interferometer, and thereby yield A = 0.75 (i.e., Ves = 0). 
More generally, the degree of imbalance for pijki can be ex- 
pressed in terms of dap, and leads to the sum uncertainty 

A< |(l-16d'). 

The pair {A, yc} thereby defines the parameter space for 



the multipartite entanglement employed in our experiment, 
with the entanglement parameters {A, Wc} serving as a non- 
local, nonlinear entanglement witness^^l Specifically, for 
a given value of yc, we numerically determine the bound 
AJ, for the minimal uncertainty possible for all states 

containing at most (Af — l)-mode entanglement and their mix- 
tures. For our quadripartite states iV = 4, we derive the re- 



(3) A (2) 



4^^^ 



A^^} for tripartite 



spective uncertainty bounds {A^ , A 

entangled, bipartite entangled and fully separable states, as 

functions of yc- A measurement of quantum statistics yc and 



the associated coherence A with A < A 



(1,2.3) 



thereby mani- 



festly confirms the presence of genuine M = 4 partite entan- 
glement, with refs.^l^li^ providing further details. 

We emphasize that our criterion of 'genuine' Af -partite en- 
tanglement based upon {A, y^} takes the most stringent form 
of non-separabilit>^" and excludes all weaker forms of multi- 
partite entanglement, including those which cannot be merely 
separated into two groups along any splittings of M (ref.^^l). 
The genuine Af -partite entangled state created from our ex- 
periment can only be represented as mixtures of pure states 
that all possess Af -partite entanglement (similar to the no- 
tion of k-producibility in multipartite spin models"' '^). Fur- 
thermore, our verification protocol is capable of unambigu- 
ously distinguishing genuine Af and (Af — l)-partite entangled 
states for any i\f < A^ by observing the uncertainty A below 
A[ ' . We take caution that our entanglement verification 
protocol cannot be applied for verifying the absence of entan- 
glement for the physical state p\^,' of infinite dimensions'^, as 

our analysis assumes a finite-dimensional truncated pr of pyp 
leading to lower bound entanglement'^. 

Figure l2]presents our results for quadripartite entanglement 
achieved for storage time tq = 0.2 ps. By appropriate settings 
of the interferometric phases for the writing and heralding pro- 
cesses in step (1), we prepare a correlated atomic state where 

(A) (A) 

the coherent component pc of p\^ is established with well- 
defined phase relations among the four ensembles {a, b, c, d} 
as for the ideal VF-state in Eq. [T] After a delay tq for step (2) 
quantum storage, the atomic state Py^ (to) is mapped to the 
photonic state p\^ (tq) in step (3) by way of read pulses. For 
step (4), the phases {/3i, /?2, /^s}^ of the verification interfer- 
ometer are optimized to match the orientations of the projec- 
tors {rij } to the stable photonic phases {(p[, 02, "^3} '^^ Pw ■ 
The photon probabilities for qij^i and Pijki are then recorded 
as in the upper and lower panels of Fig. [it. 

We first investigate off-diagonal coherence for the pur- 
portedly entangled atomic and photonic states, p\y (t) 

and p\^ (t), in Fig. ba. As the bipartite phase (32 
is varied, we observe interference fringes in the outputs 
{piooo,foioo,Pooio,Poooi} of the verification interferometer, 
and hence a variation in A, that results from the coherence 
between the bipartite entangled components of Pyyir) for 
the modes {02,62} and {c2,d2}- Furthermore, for optimal 
phase settings of (32, the observed values of A fall below the 



(3) A (2) 



4^' 



AJ,^''} (red, green, purple bands) for the 



bounds {A^-'^A 

measured value of yc and thereby certify that a fully quadri- 



a 0.75 




separable states (A < A^ ). Quantitatively, we find the un- 
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FIG. 2: Quadripartite entanglement among four atomic ensem- 
bles, a, Phase coherence between the bipartite entangled pairs of 
the quadripartite state. We observe an interference fringe between 
the bipartite modes {02,62} and {02,0(2} of the full quadripartite 
state (black points) as a function of bipartite phase fi^- As a control 
experiment, we also generate a 'crossed' quantum state p\ whose 
bipartite components for the split between {a, b} and {c, d} are in 
principle distinguishable (orange points and line), b, Exploring the 
entanglement space {A,yc} for quadripartite states. By controlling 
the spin-wave statistics for the state of the four ensembles, we ob- 
serve transitions from quadripartite, to tripartite, to bipartite entan- 
gled states, and to fully separable states (black points). Results for 

(A) 

the 'crossed' quantum state p\r are shown by the orange points and 
line, c, Expanded view of entanglement parameters near the origin 
for { A, 2/c}- We also display results for entanglement thermalization 
{ A'"^' , j/c } of the spin systems p}-, {pQ ) by the black (olive) 
dashed line. The red, green and purple bands represent the mini- 
mum uncertainties for three-mode (AJ, ) and two-mode entangle- 
ment (AJj ), and for fully separable states (A^ ). In all cases, error 
bars for the data points reflect the statistical and systematic uncer- 
tainties. The thicknesses of the {red, green, purple} bands from the 
central lines correspond to the ±l-s.d. uncertainties of the bounds 
{ A|^ , Ajj , ^f, }, due to losses and imbalances of the verification 
interferometer. 



partite entangled state p)JJ (and thereby p\^ ) has been gen- 
erated. In particular, this state verifiably possesses stronger 
quantum correlations than any admixtures of tripartite entan- 
gled (A < A^^^), bipartite entangled (A < Af ^), and fully 



(3) 



0.19] 



0.01 



certainty violation A = O.lOiO.Ol < A 

at j/c — 0.06 ± 0.02 and ^2 = 175°. This is in reasonable 

correspondence to the result Ath = 0.08 from our theoretical 

model. 

The quadripartite entanglement observed in Fig. l2h arises 
from the intrinsic indistinguishability of probability ampli- 
tudes for creation of one collective excitation \s)^ among the 
four ensembles, as implemented by the heralding measure- 
ment Ilh for fields 71 = (ai, bi,ci,di}. As a control exper- 
iment, we reconfigure the heralding interferometer such that 
the fields {ai,6i} and {ci,di} are combined with orthogo- 
nal polarizations (BSx) immediately before the heralding de- 
tector Dfi in Fig. fTJi. In this case, path-information could 
in principle be revealed up to the bipartite split of the en- 
sembles {a, b} and (c, d} by discriminating the polarization 
state of the heralding photon. Therefore, the heralding mea- 

(A) 

surement IIx prepares a 'crossed' atomic state p^ with no 
coherence shared between the ensembles {a, b} and (c, d}. 
We observe an absence of interference as (32 is varied in Fig. 
l2k (orange points). However, this modified heralding pro- 
cess fix preserves bipartite entanglement within the ensemble 
pairs {a, b} and (c, d}, which explains our observation of the 
uncertainty A = 0.533 ± 0.005 (close to the predicted re- 
suit Ajjj ' = 0.536 (orange line)) reduced below the 1-mode 



bound A 



(1) _ 



0.70 



+0.01 



for yc = 0.07 ± 0.01, thereby sup- 



porting the presence of bipartite entanglement for the ensem- 
ble pairs {a, b} and (c, d}. 

Next, we characterize p\^' (and thereby pj^ ' ) over the full 
parameter space of (A, yd- In a regime of weak excitation for 
each of the ensemble-field pairs (e, 71}, with e = {a, b, c, d} 
and 7i — (oi, &i, ci, di}, the heralded state p\y includes 
small, higher-order components and is approximately given 
by 



p\^\T^Q)c,il-30\W)AW\ + 



a(^) ;, 



3epi1^+Or), (3) 



where /5>2 includes spin-waves with two or more quanta in 
the set of four ensembles as well as uncorrected excitations 
due to atomic noise. For excitation probability ^ ^ 0, a 
heralding event at Dh leads to a state with high fidelity to 
I W^)^ in Eq. [Tlstored in the four ensembles. However, for in- 
creasing ^, higher-order terms with multiple excitation num- 
ber become important, leading to modifications of the spin- 

(A) 

wave statistics for the heralded state p\^ and thereby to the 
entanglement parameters (A, j/c}- For £, ^ 1, we find theo- 
retically that (yc)th — 8^ and Ath — 9^ (Eq. [3]), excluding any 
external noise sources and assuming perfect spatio-temporal 
overlaps for the quantum fields 72 = (02, &2, 22,^2} in the 
verification interferometer Hence, by varying the excitation 
probability ^ via the overall intensity for the write beam, we 
can simultaneously adjust the statistical contamination j/c and 
sum uncertainty A of the entangled states {py^ , PwS across 
their entanglement spaces. 

This procedure is employed in Fig. l2b to parametrically 
increase {A, yd in tandem by increasing ^. As the statistical 



contamination jjc is raised from j/c — in the quantum domain 
of one-excitation to the classical regime of multiple excita- 
tions with yc — 1, we observe the transitions of the directly 

-Si) 



mediated by the initial atom-photon {e,7i} correlations^. 
Despite the intrinsically low preparation probability ph and 
losses in the heralding channels, the resulting quadripartite 

{A) 



measured photonic PF-states jr^ (black points) from fuUy entangled state p^ stored in the four ensembles has high 
quadripartite entangled to tripartite entangled (from A < Af ^ fideUty with the ideal W^-state in Eq. [I] namely F^^' = 



to A > Ap^), to bipartite entangled (from A < A^^ to A > 
A^ ), and finally to fully separable states (from A < Aj, 

to A > A^ ). In comparison to our former work on pho- 
tonic VK-states via coherent splitting of a photon'"*, the her- 
alded atomic and photonic VK-states {py^ ,Pw} offer qual- 
itatively richer statistical passages through the entanglement 
spaces {A, yd- Here, the quantum coherence of {pI^/ , /5J^ } 
evidenced by A is inherently linked to the statistical character 
of {ij^TP^} expressed by j/c due to increasing decorrela- 
tions between the heralding fields 71 and the excitation statis- 
tics of the ensembles with increasing ^. 

(A) 

For f ^ 1, the coherent contribution pi of the delocal- 
ized single quantum strongly dominates over any other pro- 
cesses of the full quadripartite state p\^' including vacuum 
components and multiple spin-wave excitations arising from 
the uncorrected excitations of /5>2 (Eq. pL With a heralding 
probability p^ ~ 3 x 10"'' (corresponding to ^ ~ 5 x 10^^), 
we achieve the smallest entanglement parameters of A™'" = 



(A) 



[A) 



0.07: 



-0.01 



and y™" = 0.038 ± 0.006 for the generated quadri- 



partite entangled states. These entanglement parameters are 
suppressed below the closest 3-mode boundary A^ by more 
than ten standard deviations (s.d.). Furthermore, because the 
mapping of quantum states from matter to light cannot in- 
crease entanglemenrl our measurements of p^ unambigu- 
ously provides a lower bound of the quadripartite entangle- 

(A) 

ment stored m /5^ . Thereby, the observed strong violation of 
the uncertainty relations for {A™", y™"} categorically certi- 
fies for the creation of measurement-induced entanglement of 
spin-wave excitations among four quantum memories, as well 
as for the faithful and coherent transfer of the stored quadri- 
partite entangled states to an entangled state of four propagat- 
ing electromagnetic fields. We emphasize that our analysis 
makes use of the full physical state {/5J^ , /5|y } including the 
vacuum component po and higher order terms p>2, and does 
not rely upon a spurious post-diction based upon a preferred 
set of detection events. 

As shown by the black and orange curves in Fig. l2j our 
observations correspond well to a theoretical model of the 



heralding p^^ 



^^l^{^I^Ul,Ja''^ 



[/write ), readout p|^^ 



Tr^([/j^jijPj;j, JTread), and verification steps of our experiment, 

with pg = \g){jj\. In addition to Eq. bI the model includes 
the effects of atomic fluorescence, background noises, and fi- 
nite efficiencies, as well as various imperfections in the mea- 
surements. 

In terms of state fidelity, our approach for multipartite en- 
tanglement generation compares favorably to matter systems 
utilizing local interactions (e.g., trapped ions^"' ^^). Instead of 
direct interactions, the atomic M^-state in our experiment is 
generated in a heralded fashion by a quantum measurement 



1-3^ 

(-4) _ 
th ~ 



(M^aIpi^ |W^yi)- Indeed, we derive a fidelity F^ 

from Eq. bl and thereby deduce a theoretical fidelity F^ 

0.98. Experimentally, because the measurements {11^'^^ of 
the sum uncertainty project the input state p^ into the four or- 
thonormal Ty-states |M^,;)i,, we can attribute a lower bound for 
the entanglement fidehty F^-^) = {Wi\p4\Wi) = piF^^^^ 

for the heralded atomic state p\^ , with \Wi) G {|Wi)t,}. 
Specifically, the single-excitation probability, pi, and the con- 
ditional fidelity, F{ = (Wi|p![ Wi), for the atomic en- 
sembles are determined from the respective inferences of 
the spin-wave statistics (via the retrieval efficiency r/,ead and 
He), and of the coherences (via A). For the parameters 
{A™'", y™"} and ryread = 38±4%, we deduce an entanglement 



0.9 ±0.1 ~ F, 



{A) 



fidelity F(^) =pi(y'i(i-A) + i 
for the stored atomic state (Appendix). 

Apart from the creation of a novel multipartite entangled 
state of spin-waves with tunable quantum statistics, an impor- 
tant benchmark of our quantum interface is the transfer effi- 
ciency A of the quadripartite T4^-states from matter to lighpl 
Since no known measure applies to our case, we tentatively 
define the entanglement transfer A = pii) J p{A) ]^y j-jjg j-^. 
tio of the respective physical fidelities F^-^) and F'^'^^ for the 
atomic and photonic states. In particular for ^ ^ 1, the 
photonic fidelity can be approximated by F^^' ~ Jyread-Flh ' 
which thereby gives Ath — ?yi-ead dictated by the retrieval 
efficiency reread- For the minimal entanglement parameters 
{A™",y™"} in Fig. ^, we find F^'') = 0.35 ± 0.02 and 
obtain an entanglement transfer A = 0.39 ± 0.04, similar to 
the observed Thread- While fidelity is an often used measure, 
we emphasize that F^''^' cannot be used to set a threshold for 
entanglement for p\^ , since p\y can exhibit multipartite en- 
tanglement for any F^'^> > 0. 

To investigate the dynamical behavior of the observed 
quadripartite entangled states, we study the temporal evolu- 
tion of multipartite entanglement stored in the atomic ensem- 
bles as a function of time r following the heralding event 
at r = 0. Decoherence for the atomic W^-state is gov- 
erned by independent dephasing processes in each ensem- 
ble, in which the imprinted atomic phases in the collective 
state jsj) evolve independently due to atomic thermal mo- 
tion, which transforms the initial collective atomic state into a 
subradiant state^** that becomes increasingly uncorrected with 
the heralding fields 71 — {oi, 61, ci, di}, for each ensemble 
e = {a, b, c, d}. The subradiant state contributes to a reduc- 

(A) (A) 

tion in the coherent component pc of p^ as well as to a 

build-up of uncorrected atomic noise pS^j relative to pc ' . 
The net effect is an increase of both entanglement parameters 
{A, He} with a typical time-scale of r,,,, = l/(|(5fc|W(i) ~ 17 
ps, where Sk is the momentum transfer to the spin-waves 
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FIG. 3: Temporal decay of coherences stored in four atomic ensembles, a, Evolution of the photon probabilities {piooo,Poioo,Pooio,Poooi} 

for occupying the output modes of the verification interferometer versus storage time r. For readability, the heights of the bars are shown in 
accord to the color convention of inset i. Error bars, shown as grey squares, reflect the statistical uncertainties for each point, b, Photon 
probabilities {piooo,Poioo,Pooio,Poooi} from our theoretical model, which assumes a memory time determined from the temperature of the 
cold atomic samples and the net momentum transfer to the atomic spin-waves (Appendix). 



and Vd is the average velocity of the atoms^^. Eventually, the 
growth in {A(r), ydT)} leads to time-dependent losses of en- 
tanglement, marked by successive crossings of the boundaries 

setby{Af\Ai'\AW}. 

Fig. [Slillustrates the temporal reduction in the overall co- 
herence d of the full quadripartite state. Operationally, the 
loss of coherence is observed in terms of decrease in imbal- 
ances among {piooo,Poioo,Pooio,Poooi} as a function of stor- 
age time r, and hence to an increase in A. The behavior of 
the experimentally observed photon probabilities in Fig. l3t 
results from the progressive decay of the initial coherence for 
p\^\Ti) atT, = 0.2 ^s for which V;ff(T,) = 4d = 0.95±0.02, 

evolving then to V^iJ = O.lOta^o for the final state /)|^^(r/) 
measured at ry = 36.2 /is. The observed evolution is in good 
agreement with our theoretical model of the photon probabil- 
ities shown in Fig. [3]3. The spin-wave statistics are similarly 
modified by phase decoherence leading to an increase of -ijc, 
from ydTi) = 0.03 ± 0.01 to ydTf) = 0.74 ± 0.34. 

Finally, in Fig. |4]we examine the dissipative dynamics of 
entanglement for the quantum memories of four ensembles 
via the evolution of both {A, yd- We observe the passage of 



the initial quadripartite entangled state p^ (r^) through vari- 
ous domains, evolving from M-partite entanglement to (A/ — 
1) -partite entanglement at memory times t = t„, . The 
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crossings of the bounds { A^ , Aj^ , Aj^ } occur at r, 

lis, Tm =21 lis, and Tm — 24 iis, respectively. In addition, 
we find that the measured entanglement parameters evolve in 
accord to the simulated dynamics derived for p^^ (t) from 
our theoretical model (solid line). In Fig. ^ we show the 
projection of Fig. pb into the plane of {A, y^}, which directly 
reveals the losses of entanglement, albeit without any tempo- 
ral information. 

An interesting extension of our work is to relate the charac- 
terization of multipartite entanglement by way of {A, yd as 
in Fig. |2]to the thermal relax ation of entanglement for quan- 
tum many-body systema^^E] We consider a ferromagnetic 
(J > 0) Heisenberg-like model i?^ = Hu + Hp and a Lipkin- 
Meshkov-Glick model i^LMC = -{ E(»j) (^^'^ ' si'^ + ^i"^ ' 



:0■)^ 
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(J) Mi) A(i)-| 



Sy) + h, J2i Si', for four spins {S^'> , S^y> , S^"> } on a tetra- 
hedron. Here, if^ includes the standard Heisenberg Hamil- 
tonian H^ = --{ J2{^,J) S^"^ ' S'-''^ + ^z E^ Sz\ in addition 
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FIG. 4: Dissipative dynamics of atomic entanglement, a, Time-evolution of the entanglement parameters {A(r), j/c(t)} versus r for the 
multipartite quantum state. We observe the crossing of the boundaries defining the minimum uncertainties for 3-mode (red surface, AJ^ ), 



(l)^ 



2-mode (green surface, A^ ) entangled states, and 1-mode (purple surface, Aj^ ) separable states. For clarity, the data points and theoretical 
curve are colored to indicate the various entanglement orders for the quadripartite (black), tripartite (red), bipartite entangled (green) states, 
and fully separable states (purple). The projections of the data points into the planes {yc, t) and (A, r) are shown as gray points to display 
the individual passages of the quantum statistics and coherences, respectively, b, Projection of entanglement dynamics onto the (A, yc) plane. 
The various entanglement transitions from a are shown for our measurement and theory. Error bars for the data represent 1/e-errors from a 
conservative optimization analysis for evaluating {A, j/c} that reflects the parameters' statistical and systematic uncertainties. The curves are 
from a theoretical model that includes motional dephasing. 



to a spin-projection term Hp = 2hz\S,—S){S,—S\, which 
shifts the energy of the collective spin state jS*, — 5). For- 
mally, the entangled state p\^ in Eq. BI closely resembles the 
thermal equilibrium state pf^ ip'^^^^) of H!^ (Hlmg)- The 
eigenstates of _ff^ and TJlmg ^e the total angular momen- 
tum states defined by collective spin algebra, whose ground 
state at zero temperature (T — ^ 0) is an entangled M^-state 
similar to Eq. [Tlfor magnetic field h^ — J/2. By redefin- 
ing {| I), I t)} i-^ {|0), |1)}, we calculate the Gibbs state 

p(H) = (l/Z)e-^H/fcBT (p(LMG) ^ (l/Z)e-f^LMG/fcBT) ^jjij 

partition function Z as a function of {A^'^\yi '} derived 
for this system fl^ (flLMc)- The results are shown in Fig. 
|2j; by the black (olive) dashed lines for p(f' (p^^^^)). The 

(A) 

Statistical character of pyy' (black points) for our system of 
four ensembles follows closely the thermalization process rel- 
evant for the low-lying states of _ff^ (black dashed line) for 

Vc ^ 0.2, whereby the higher-order spin-waves p^J are pop- 
ulated by thermal excitations. While intriguing, such an anal- 
ogy is incomplete for our system, because pyy' is comprised 
of four NA,e dimensional quasi-bosonic excitations instead of 
the four single-spins {| |), | t)} of p^^K p'^^^'^K Neverthe- 
less, our methodology for multipartite entanglement charac- 
terization is clearly applicable to quantum spin systems and 
may help to advance our understanding of the link between 
off-diagonal long-range or der an d multipartite entangled spin- 
waves in quantum magnets^^ED. 

In conclusion, our measurements explicitly demonstrate a 
coherent matter-light quantum interface for multipartite en- 
tanglement by wayof the operational metric of quantum un- 
certainty relationg^^GD High-fidelity entangled spin-waves 



are generated in four spatially separated atomic ensembles 
and subsequently coherently transferred to quadripartite en- 
tangled beams of light. The quantum memories provided by 
the four atomic ensembles are individually addressable and 
can be readily read-out at different times for conditional con- 
trol of entanglement distribution and connection"*. By way 
of recent advances from other groups (Appendix), the short 
memory times for our work could be significantly improved 
beyond one second. In addition, measurement-induced gen- 
eration of entanglement could be augmented by the reversible 
mapping of entangled photonic states to matter". 

Other possibilities include the creation of yet larger mul- 
tipartite entangled states with efficient scaling^. For exam- 
ple, quadripartite entangled states of ensemble sets {a, b, c, d} 
and {a' , b' ,c' , d'} could be extended by swapping between en- 
sembles {a, a'} to prepare a hexapartite (6-partite) entangled 
state for ensembles {b, b' , c, c', d, d'} (see Appendix). Gener- 
alization of such entanglement swapping processes may help 
to prepare a single macroscopic entangled state which fully 
occupies an entire network, for exploratory goals of observ- 
ing the extreme non-locality of M^-states^" and entanglement 
percolation^ . In addition, the expansion of multipartite quan- 
tum networks with increasing N offers a promising bottom- 
up approach for the study of phase transitions in strongly- 
correlated systems'^ '^. Finally, the entangled spin-waves 
comprising the network can be used to implement a quantum- 
enhanced parameter estimation scheme for detecting an un- 
known phase shift of tt for one component of p\^. with effi- 
ciency beyond any separable state, which could find applica- 
tion to quantum metrology. 
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Appendix 

Experimental details. The experiment consists of a 22 ms 
preparation stage and a 3 ms period for operating the quantum 
interface in Fig. [Tlwith a repetition rate of Ri = 40 Hz. For 
the preparation, we load and laser-cool Caesium atoms (peak 
optical depth « 30) in a magneto-optical trap for 18 ms, af- 
ter which the atoms are further cooled by optical molasses 
(Td ~ 150 /.iK) and optically pumped to \g) for 4 ms. During 
this time, a phase-reference laser (F — 3 ^ F' ^ 4' tran- 
sition) also propagates through the atomic ensembles for the 
active stabilization of the verification interferometer in Fig. [T]; 
via ex-situ phase-modulation spectroscopy '"*, which does not 
affect the operation of the quantum interface. Concurrently, 
dense Caesium atoms in paraffin coated vapor cells located 
at the heralding and verification ports are prepared to the re- 
spective ground states \g) (\s)) for filtering the classical lasers 
scattered into the quantum fields 71 (72). 
Quantum interface. For the quantum interface to func- 
tion during the 3 ms window, in step (1) 20-ns wide writ- 
ing (red-detuned 5 = 10 MHz from \g) — |e) transition) and 
reading pulses (resonant with \s) — |e)) are applied sequen- 
tially to the ensembles e, synchronized to a clock running 
at i?c — 2 MHz. This process creates pairwise correlated 
excitations^ between the collective atomic modes |s)e of the 
ensembles e and the optical fields 71 (5 ~ 10 MHz below 
\s) — |e)). Photodetection of a single photon for the com- 
bined fields 71 at the output of the heralding interferometer 
effectively erases the which-path information for 71 , and im- 
prints the entangled spin-wave p\^' (Eq. bI onto the ensem- 
bles {a,6,c, d} via Txh{t\hUlMj^g^^Uv,ni<,)- The heralding 
event at Dh triggers control logic in Fig. [Tk which deactivates 
intensity modulators of the write (IMwrite) and read lasers 

(A) 

(IMiead) for the quantum storage of p\^' in step (2). After 
a user-controlled delay t, step (3) is initiated with strong read 
pulses (Rabi frequency 24 MHz) that address the ensembles in 
Fig. fTt and coherently transfer the entangled atomic compo- 
nents {a, 6, c, d\ of p\^ (t) one-by-one to propagating beams 
72 — {a2,b2,C2,d2}, comprising the entangled photonic 
state p'^\t), via the operation p'^^ = TvAiUl.,ip\^^U,s^i). 
Here, Tta traces over the atomic systems which are later 
shelved into the ground states |^). The retrieval efficiency 
?7read is collectively enhanced for large Na (ref.l^, leading to 
jyread — 0-38 ± 0.06 in our experiment. 
Spin-wave quantum memories. The quantum information of 
the entangled state for Eq. [Tlis encoded in the quantum num- 
bers of spin-waves (collective excitations) for the pseudo-spin 
of the hyperfine ground electronic levels 65*1/2! F = 3,F ~ 4 
in atomic Caesium. The fluorescence images shown in the in- 
set of Fig. [Tk arise from excitation by the writing and read- 
ing beams with 1 mm separations and 200 p,m waists. The 
geometry of actual collective excitations for the four ensem- 
bles {a, b, c, d} is defined by the spread functions of the imag- 
ing systems for the fields 71, 72, which have waists of 60 /im 
for each ensemble consisting of a cold cloud of Na.c ~ 10® 
Caesium atoms. We use an off-axial configuration-' ' for ad- 
dressing each ensemble e with an angle 6 = 2.5° between 
the classical and nonclassical beams, that creates spin-waves 



|s)e associated with wave-numbers 6k — /cwrite — ^71 for each 
e. These spin-waves are analogous to other types of collec- 
tive excitations in many-body systems, such as magnons and 
plasmons, and the spin-waves can be converted to dark-state 
polaritons for the coherent transfer C/iead of entanglement. For 
the phase-matching configuration and temperature of our en- 
sembles, the memory times {rm ,t„i ,t„i } in Fig. Ware 
dominantly determined by the motional dephasing of the spin- 
waves |s)(p2l With thermal velocity of w^ ~ 14 cm/s, we the- 
oretically determine a memory time t„i = ^^"^'^^f^,™-),, =17 
fis. On the other hand, the ground-state dephasing due to in- 
homogeneous broadening is expected to be > 50 us in our 
experiment, inferred from two-photon Raman spectroscopy. 
Derivation of entanglement fidelity. We derive here the ex- 
pression for the lower bound entanglement fidelity f('^) = 

" (c) 

PiFi. We start by noting that the projective measurement 11^ 
for A gives the conditional fidehty Fi of the input Pr onto one 
of four orthonormal M^-states, |Wi)„ = |Wi)u, for example 



1 1000) 



=»/?. 



1 10100) + e'/^^ (10010) + e'^^ |0001)). Hence, we 



can define A = 1 — Ff — J2i=2 ^i ^^ terms of the respective 
overlaps Fi. Because of the orthonormality X]i=i ^i — 1' the 



sum uncertainty is bounded by A > 1 — Ff — (1 — FiY, 
whereby we obtain Fi > ^/^(i — A) + ^. Finally, by 



combining the probability pi for exciting one spin-wave dis- 
tributed among the four ensembles, we access the lower bound 

fidehty F^^) > pi ( J^{^ -A) + i ) for the heralded atomic 

(A) 

State p\y ■ In principle, the imbalances in the interferome- 
ter rotate the projectors into non-orthonormal sets'^. How- 
ever, the measured losses and the beam-splitter ratios are suf- 
ficiently balanced such that any changes in F(^) due to mod- 
ified projectors are well within the uncertainties of the data. 
Prospects for improving memory time and matter-light 
transfer efficiency. By operating the clock speed at Re -^ 10 

(3) 
MHz and Tm « 20ps, we can prepare hexapartite (M = 6) 

entanglement with probability of 3z?7readPft,/8 ~ 10^^ by 

(A) 

connecting two quadripartite states p\^ , where the enhance- 
ment factor z is 400 (ref.-'^). However, the memory times 
{tL ,Tm ,Tm } in Fig. Wand the entanglement transfer A 
from matter to light limit our capability to scale the mul- 
tipartite entanglement beyond M > 6 by w ay of quan- 
tum control and swapping of entanglemenP^EH with our cur- 
rent experimental parameters. The prerequisite storage tech- 
niques for suppressing both the internal and motional spin- 
wave dephasings can be extended for t,„ with advances in 
ensemble-based quantum memories^'* ■'^. Recent experiments 
with single ensembles have achieved coherence t imes up to 
Tm — 1-5 seconds in quantum degenerate gases^^El albeit 
with efficiencies < 1%. The transfer efficiency can also be 
increased to Ath — 0.9 by enclosing the ensembles with high 
finesse cavities^. System integrati ons b y way of atom-chip 
technology and waveguide coupling^"^' hold great potential 
for scalability given the strong cooperativity and the long 
coherence'*'. In addition, our experiment opens the door for 
future theoretical studies of complex repeater architectures for 
multipartite systems, beyond traditional one-to-one networks. 
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